AD-A182  388 
UNCLASSIFIED 


FURTHER  RESULTS  IN  BEND-BUCKLINQ  ANAL  VS  IS  OF  RING  1/1 

STIFFENED  CYLINDRICAL  SHELLS(U)  NAVAL  SURFACE  HEAPONS 
CENTER  SILVER  SPRING  HD  H  HOUSSOUROS  AUG  88 
NSMC/TR-8S-226  F/G  28/11  NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDSJ963-A 


AD-A182  380 


*•  <  «t*  *  i». 


NSWC  TR  86-326 


OJK  FILE  CD 


FURTHER  RESULTS  IN  BEND-BUCKLING  ANALYSIS 
OF  RING  STIFFENED  CYLINDRICAL  SHELLS 


BY  MINOS  MOUSSOUROS 


RESEARCH  AND  TECHNOLOGY  DEPARTMENT 


AUGUST  1986 


Approved  for  public  release;  distribution  is  unlimited. 


i/'Orv  or  f 


DT1C 

^ELECTE 
K  JUL1  61987 


NAVAL  SURFACE  WEAPONS  CENTER 

Dahlgren,  Virginia  22448-5000  •  Silver  Spring,  Maryland  20903-5000 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  7H iS  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a  SECURITY  CLASSIFICATION  AUTHORITY 


AtPiifo 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  MARKINGS 


2b  DECLASSIFICATION  '  DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

NSWC  TR  86-326 


3  DISTRIBUTION  .'AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
is  unlimited. 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a  NAME  OF  PERFORMING  ORGANIZATION  6b  OFF.CE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 

(If  applicable) 

Naval  Surface  Weapons  Center  Code  R14 


6c  ADDRESS  (City,  State,  and  ZIP  Code) 

10901  New  Hampshire  Avenue 
Silver  Spring,  MD  20903-5000 


8a  NAME  OF  FUNDING  ,  SPONSORING 
ORGANIZATION 


7b  ADDRESS  (City.  State,  and  ZIP  Code) 


8b  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 


8c  ADDRESS  (City,  State,  and  ZIP  Code) 


1 1  Title  (Include  Security  Classification ) 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO. 

NO 

NO 

ACCESSION  NO 

62633N 

F33327 

SF33327 

R19BA 

Further  Results  in  Bend-Buckling  Analysis  of  BirfgStif f ened  Cylindrical  Shells 


12  PERSONAL  AUTHOR(S) 

Moussouros,  Minos 


13a  type  Of  report 


16  SUPPLEMENTARY  NOTA’ION 


14  DATE  OF  REPORT  (Year,  Month.  Day)  15  PAC 
1986,  August  5 


COSATI  CODES 


13 


J  Sr, 


18  9B8JECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

group _ Subgroup  I  bendibuckling  ^  finite  element  i  J 

13  I  cylindrical  shells  j 


’9  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number ) 

This  report  presents  finite  element  bend-buckling  analysis  results,  including  short 
wavelength  (local)  buckling,  and  an  analysis  of  the  surface  of  a  nearly  circular  cylinder 
through  a  double  Fourier  series.  /T 


i  - "v  ^  !  - - -  — 


20  D>STR'8UT'ON  AVAiLAB'UTy  OF  ABS’RAL’  2'  ABSTRACT  SECURITY  ClASSiF'CANON 

□  UNCLASSIFIED  UNL  M ’ED  QD  SAME  AS  RPT  □  DTiC  USERS  UNCLASSIFIED 


22a  NAVE  OF  RESPONS  0.E  'ND-J'Dua,  22b  TELEPHONE  (Include  Area  Code)  22 c  OFFICE  SYMBOL 

Minos  Moussouros  (202)394-1681  Code  R14 


00  FORM  1473,  84  VAP  83  APR  ed  T  on  T*,  oe  »seo  until  exhausted  SECURITY  CLASSIFICATION  OF  th.s  RAGE 

a  i  itt'cr  e<j. font  are  3Dso'ete 


ft  U  S.  CoMftlMM  Pr.bt.n,  Office  407  0*4 

UNCLASSIFIED 


NSWC  TR  86-326 


FOREWORD 


This  report  presents  finite  element  bend-buckling  analysis  results, 
including  short  wavelength  (local)  buckling  and  an  analysis  of  the  surface  of  a 
"nearly  circular  cylinder”  through  a  double  Fourier  series. 

This  study  used  the  nonlinear  finite  element  program  ABAQUS.  It  was 
performed  under  an  ongoing  program  aimed  at  realistic  modeling  of  shell 
structures  and  capturing  all  the  developing  phenomena  during  static  buckling  of 
circular  cylindrical  shells  subject  to  either  bending  moments  or  axial 
compression. 

This  work  was  supported  by  the  Office  of  Naval  Technology  through  the  Naval 
Surface  Weapons  Center  Block  Program,  "Undersea  Weapons,  Warheads  and  Fuzes" 
(D.E.  Phillips). 


Approved  by: 

- 

KURT  F.  MUELLER,  Head 
Energetic  Materials  Division 
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INTRODUCTION 


The  ability  to  analyze  cylindrical  shells  subjected  to  bending  loads  was 
the  topic  of  an  earlier  study. 1  In  that  report,  three  unstiffened  long 
circular  cylindrical  shells  were  analyzed  by  the  nonlinear  finite  element 
program  ABAQUS.  The  computational  results  were  then  compared  with  experimental 
data.  The  agreement  was  good.  This  report  presents: 


1.  Further  analysis  results  of  circularly  perfect  models  subject  to 
end-bending  moments.  This  is  a  follow-up  of  References  1  and  2,  including  short 
wavelength  buckling  for  ring-stiffened  cylinders. 


2.  An  analysis  of  the  surface  of  a  "nearly"  circular  cylinder,  given  its 
Fourier  components;  or  inversely,  construction  of  the  Fourier  coefficients  from 
given  data  representing  the  surface. 23,24  ye  intend  to  be  able  to  generate 
imperfect  geometries  from  a  perfect  model  and  easily  generate  the  topological 
model  employed  in  the  ABAQUS  general  nonlinear  finite  element  computer 
nmomm  25 


program. 


Furthermore,  this  report  extends  the  list  of  references  in  References  1  and 
2.  References  3  through  13  represent  additional  works  on  bend-buckling  analysis 
of  cylindrical  shells.  References  14  through  22  are  additional  general  works  on 
analysis  of  shells.  Reference  23  deals  with  the  fast  Fourier  transform 
algorithm  which  was  not  used  here. 


The  methods  applied  here  in  connection  with  the  surface  representation^ >24 
of  a  "nearly"  cylindrical  shell  are  described  in  Reference  24.  Appendix  A 
contains  the  mathematical  details  of  the  method  employed  in  this  report. 


Appendix  B  contains  the  actual  FORTRAN  listing  of  the  subroutine  to  be  used 
with  ABAQUS  in  these  computations. 


WWW 


rag* 
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METHOD 


FINITE  ELEMENT  MODELING 

The  analysis  of  a  circular  or  near-circular  cylindrical  shell  subject  to 
end  bending  moments  will  be  addressed  here  numerically.  For  this  purpose  the 
nonlinear  finite  element  program  ABAQUS^  is  used. 

In  this  work,  we  establish  first  our  global  Cartesian  frame  of  reference 
(x,y,z,  Figure  1).  The  shell  axial  direction  is  along  z.  The  vertical 
direction  is  Indicated  by  x,  while  the  transverse  y  is  along  an  axis  pointing 
out  of  the  paper  towards  the  reader.  Again  only  180°  of  the  structure 
peripherally  and  L/2  axially  are  discretized  with  finite  elements.  If  the 
cylinder  is  not  exactly  circular,  then  Fourier  analysis  of  the  kind  indicated  in 
Appendix  A  will  produce  the  various  dominant  harmonics  axially  and 
peripherally.  In  such  cases,  it  may  be  necessary  to  analyze  360°  of  the  shell 
peripherally  and  even  its  entire  length.  Here  it  is  assumed  that  both  cylinder 
and  loads  are  symmetrical  for  this  simplification  to  be  correct. 

We  discretize  the  cylindrical  surface  by  ABAQUS  S8R  shell  elements  with  3 
integration  points  across  the  thickness.  This  means  that  both  outside,  inside, 
and  midsurfaces  of  the  shell  are  considered,  when  stresses  are  computed  later. 
The  element  consists  of  8  nodes,  4  corner  and  4  midside  ones,  and  employs 
reduced  integration  to  calculate  the  stiffness  matrix.  The  program  allows  for 
geometric  and  material  nonlinearities.  Furthermore,  the  constitutive  equations 
consist  of  a  von  Mises  yield  surface  with  the  associated  flow  rule.  The 
hardening  rule  considered  here  is  isotropic. 

Peripherally,  the  discretization  allows  25  nodes  (including  midside  ones) 
over  180°.  Axially,  depending  on  the  length  analyzed,  we  have  a  variety  of 
total  nodes. 
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BOUNDARY  CONDITIONS  AND  APPLIED  LOADING 

The  boundary  conditions  are  given  in  the  global  system.  (See  Figure  1.*) 
They  are: 

1.  Symmetry  about  the  x-z  plane. 

At  8  =  0  and  0  =  tt  we  must  have: 
v  =  V’jc  =  v?z  =  0 

2.  Symmetry  about  the  x-y  plane  at  midlength  (x  =  L/2). 

W  =  V’x  =  <fiy  =  0 

3.  We  must  also  remove  the  vertical  rigid  body  motion  (x  axis)  and 
apply  our  external  loading  in  the  form  of  a  rotation  about  the 
global  transverse  axis  y.  These  two  aspects  have  been  explained  in 
length  in  References  1,  2,  and  25,  and  are  omitted  here. 


APPLICATIONS  AND  RESULTS 

This  report  contains  results  on  two  types  of  cylindrical  shells, 
unstiffened  and  ring-stiffened.  The  unstiffened  shells  represent  variations  of 
the  Jisra  models  reported  in.^>^  The  ring-stiffened  models  represent  results 
of  shells  already  tested  in  purely  axial  compression^  and  can  only  be  used 
for  further  study,  since  no  other  ring-stiffened  data  suitable  for  marine 
applications  are  available. 

In  this  report  we  analyzed  again  models  10A,  16A,  and  20A  of  Jisra^>2  and 
they  will  be  referred  to  as  models  10AN,  16AN,  and  20AN,  respectively.  They 
differ  from  models  10A,  16A,  and  20A  because,  in  this  case,  the  analyzed  span 
Lr  is  shorter,  as  indicated  in  Figure  2.  Tables  1  and  2  give  the  model 
characteristics . 

Figure  3  is  a  comparison  of  the  analysis  of  models  10A  and  10AN.^»^ 
Furthermore,  the  included  tabular  form  can  facilitate  the  comparison.  Figures  4 
and  5  report  results  for  models  16AN  and  20AN  compared  with  previous 
computations.  The  disagreement  is  more  pronounced  for  model  20AN  than  the  first 
two.  Looking  at  Figures  4  and  5  we  also  notice  that  the  slopes  are  not  exactly 
the  same.  However,  if  more  initial  points  were  saved  at  the  beginning  of  the 
computations,  then  the  slopes  would  have  agreed.  This  is  known  from  model  10AN , 
for  which  the  initial  analyses  were  stored  for  this  purpose.  Our  first 
conclusion,  therefore,  is  that  the  thicker  and  shorter  the  cylinder  is,  the  more 
important  the  length  becomes  (reduction  in  computed  bending  moment  was  about  9%). 


*  v  =  Translation  along  global  y  axis  (laterally) 

w  =  Translation  along  global  z  axis  (longitudinally) 
<PX  =  Rotation  about  global  x  axis 

'f y  =  Rotation  about  global  y  axis 

¥>z  =  Rotation  about  global  z  axis 
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Figures  6,  7,  and  8  pertain  to  model  10AN.  In  Figure  6  the  moment-rotation 
relation  is  given.  The  table  in  Figure  6  gives  the  corresponding  total  energies 
as  deformation  progresses.  Figure  7  shows  how  the  bending  moment  is 
independently  computed  for  various  angles  of  rotation  (indicated  as  steps*). 
Axial  and  hoop  stresses  are  also  given  for  the  first  stored  step  of  the 
analysis.  Figure  8  gives  the  stresses  for  the  last  stored  step. 

Figure  9  pertains  to  model  16AN  and  is  similar  to  Figure  6.  It  is  the 
moment-rotation  relation  that  is  displayed.  Figure  10  contains  the  same 
information  as  in  Figures  6  and  9,  but  it  applies  to  model  20AN.  Figure  11  is  a 
tabular  form  of  the  axial  and  hoop  stresses  at  midlength  for  the  very  first 
increment,  i.e.  for  an  angle  of  0.948961  degs.  The  previous  given  figures 
(6-11)  only  add  to  a  better  understanding  of  the  type  of  results  to  expect  with 
the  method  so  far. 

Table  3  gives  the  ordinates  of  the  imperfect  model  20A  and  Table  4,  the 
corresponding  double  Fourier  coefficients.  We  observe  that  the  imperfections  of 
this  model  are  over  360°.  Figures  12  and  13  display  some  of  the  harmonics.  For 
example,  Figure  12  gives  the  magnitude  of  the  fundamental,  first,  and  second 
axial  harmonics  of  the  cosine-cosine  coefficients.  Similarly,  Figure  13 
displays  the  magnitude  of  the  fundamental,  first,  and  second  peripheral 
harmonics  of  the  cosine-cosine  coefficients.  The  purpose  of  this  analysis  is  to 
establish  the  software  for  quickly  producing  the  dominant  harmonics  of  imperfect 
cylindrical  shells.  It  is  known  that  preexisting  mode  patterns  will  affect  the 
final  pattern  whether  there  will  be  a  change  of  mode  or  not.  In  fact,  there 
will  be  an  amplification  along  the  "preferred  imperfection."  Such  behavior  will 
induce  buckling  and  lead  to  faster  collapse.  The  application  is  certainly 
obvious,  not  only  for  bend-buckling  but  for  other  forms  of  structural 
instabilities,  static  or  dynamic.  As  a  result,  a  FORTRAN  computer  program  was 
written  for  VAX  11/785,  implementing  Appendix  A  in  such  a  way  that  the  user  can 
either  obtain  the  Fourier  coefficients  of  given  imperfection  data  or,  given  the 
imperfection  Fourier  coefficients,  obtain  the  configuration.  Here  we  observe 
again  that  some  of  the  boundary  conditions  will  have  to  change  to  correctly 
account  for  360°  models,  and  bend-buckling  results  with  imperfections  will  be 
presented  elsewhere. 

Tables  5  and  6  display  some  of  the  geometric  characteristics  of  the 
ring-stiffened  shells  which  form  the  second  topic  of  this  report.  The  models 
are  referred  to  by  names  IC1  and  IC3.2&>27  They  were  chosen  because,  at 
least,  axial  compression  tests  had  already  been  carried  out  and  reported 
elsewhere. 26  Figure  14  represents  the  deformed  profile  of  model  ICl  at  step 
3,  increment  30,  i.e.  for  an  angle  of  rotation  of  0.231937  degs,  as  can  be 
deduced  from  Figure  16.  Figure  15  displays  the  longitudinal  and  hoop  stresses 
for  the  very  first  stored  step  (step  1,  increment  5).  We  observe  that  the  hoop 


*In  these  analyses,  the  rotation  is  increased  gradually.  During  such  a  process 
or  "step"  or  "increment,"  various  quantities  are  computed  and  stored.  For 
example,  in  Figure  6,  what  the  ABAQUS  computer  program  refers  to  as  "step  1, 
increment  15"  was  the  third  "stored  increment"  or  "step." 

The  step  size  is  not  constant  but  usually  changes  during  the  computation.  It 
increases  linearly  and  in  the  present  context  it  represents  change  of  angle  of 
rotation.  The  reported  values  represent  the  stored  values  and  not  the  actual 
number  of  steps  carried  out  by  ABAQUS. 
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stress  distribution  at  midlength  (because  this  is  the  location  for  which  it  has 
been  plotted  together  with  the  longitudinal  stresses  at  the  midsurface)  differs 
markedly  from  the  stress  distribution  of  unstiffened  cylinders  in  bending  (see 
Figure  7  for  example). 

Figure  16  relates  to  energies,  bending  moments,  angles  of  rotation,  while 
Figure  17  gives  results  of  bending  moments  from  consistent  loads.  Figure  18  is 
additional  information  for  model  IC1. 

Figure  19  plots  cross-section  flattening*  versus  mean  curvature. 

Figure  20  displays  external  work  versus  nondimens ional  bending  moment 

Figure  21  displays  cross-section  flattening  versus  M/Mcr«  Figure  22 
gives  the  distribution  of  axial  and  hoop  stresses  at  the  last  step  (step  3, 
increment  30)  graphically  and  numerically.  Finally,  Figure  23  contains  3 
plots.  The  last  two  are  moment  curvature  plots^  while  the  first  represents  the 
plot  of  the  difference  of  radial  displacements*  from  the  left  end  (OX)  of  the 
cylinder,  where  the  moment  (rotation  actually)  is  applied  up  to  midlength 
(50%).  This  is  for  the  compression  side  at  midthickness.  Observe  that  at  the 
left  end  (0%)  this  difference  is  nonzero,  because  the  radial  displacement  w  is 
not  restricted  there.  Essentially,  similar  results  (Figure  24)  are  obtained 
when  the  displacement  w  is  restricted  (w  =  o)  at  the  left  end  cross-section, 
where  the  bending  moment  is  being  applied.  This  is  the  case  in  this  situation 
because  of  end  plates.  Notice  that  in  this  case.  Figure  25  compares  with  the 
top  plot  of  Figure  23,  except  that  at  the  left  end  (0%)  the  nondimens ional 
difference  of  radial  displacements  is  zero. 

By  Figure  16  or  18,  when  the  bending  moment  is  about  21,574  K-in  (at  center 
and  end  sections  of  cylindrical  shell)  the  finite  element  program  (ABAQUS) 
predicts  short  wavelength  (local)  buckling.  Table  6  indicates  that  the  critical 
theoretical  bifurcation  bending  moment  of  the  end  panel  there  reaches  21,495.3 
K-in  with  a  theoretical  end  panel  critical  stress  of  69.5  ksi.  The  approximate 
compressive  buckling  stress  by  ABAQUS  was  about  49.18  ksi.  (See  Figure  26.)  In 
all  of  these  finite  element  calculations,  the  bending  moment  has  been  obtained 
in  two  different  ways,  i.e.  either  as: 

(1)  A  direct  output  from  ABAQUS  (See  bottom  of  Figure  16,  where  bending 
moments  are  given  in  terms  of  angle  of  rotation.)  In  the  case  of  model  IC1,  the 
critical  bending  moment  was  calculated  this  way  as  21,574  K-in. 


*Cross-section  flattening  (Aw/h)  is  defined  as: 

Aw/h  =  [  w(L/2,0)  +  w(L/2,tt)  ]/h 

More  specifically,  consider  the  radial  displacement  w  at  a  load  level  (in 
the  present  instant  rotation  level)  at  which  instability  is  assumed  to  occur 
(reference  state).  Then  following  Almroth's  suggestion,  we  construct  the 
nondimens ional  difference  of  radial  displacement  distributions  £w/h  at  the 
compression  side  for  load  Levels  higher  than  the  reference  one.  These  are 
plotted  axially  and  are  referred  to  as  "critical  parameter  for  short  axial 
wavelength  buckling."  Figures  19,  20,  21,  and  23  indicate  discontinuities  at 
the  point  of  buckling. 
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(2)  Indirectly,  based  on  the  consistent  reaction  forces  and  their 
respective  arms  from  the  position  of  neutral  axis.  (See  the  last  four  pages  of 
Figure  17.)  In  the  case  of  model  IC1,  such  an  analysis  gave  2  X  10,783.8  •  21.567 
K-in.  (Refer  also  to  Figure  18  for  such  a  summary.) 

Reference  26  (Table  5)  gives  a  much  lower  "failure  stress"  (25.6  ksi)  for  the 
center  panel,  for  the  case  of  pure  compression.  Figure  26,  used  in  conjunction 
with  Figure  16,  gives  the  axial  and  hoop  stress  distributions  for  model  ICl  at 
step  2,  increment  20,  i.e.  around  the  onset  of  buckling. 

Figure  27  displays  the  deformed  and  undeformed  profile  of  model  IC3. 

Figures  28  and  29  are  plots  of  bending  moment  versus  angle  of  rotation  and  a 
nondimens ional  curvature,  while  Figure  30  gives  the  critical  parameter  for  short 
wavelength  (local)  buckling. 

Figures  31  and  32  present  the  axial  and  hoop  stress  distributions  for  model 
IC3  for  two  different  loading  states.  Figure  31  corresponds  to  step  2, 
increment  25,  while  Figure  32  to  step  2,  increment  50.  This  means  that  Figure  31 
corresponds  to  the  onset  of  buckling  (giving  an  approximate  maximum  compressive 
stress  of  51.72  ksi)  while  Figure  32,  to  the  last  recorded  load  step  with  a 
maximum  compressive  stress  of  52.2  ksi.  Therefore,  for  model  IC3  ABAQUS 
directly  predicted  a  peak  bending  moment  of  8,689.45  K-in  (Figure  33)  or  by 
indirect  computation  (as  explained  for  model  ICl)  8,690.34  K-in  (Figure  34). 

This  can  be  seen  from  Figure  33  which  relates  energies,  moments,  angle  of 
rotation,  step,  and  increment  for  model  IC3.  Here  some  further  background  is 
necessary  to  bridge  these  results  with  the  few  experimental  data  that  are 
available.  In  1961  (see  Reference  59  of  Reference  1),  Seide  and  Ueingarten 
performed  a  bifurcation  analysis  for  a  circular  cylindrical  shell  subject  to  end 
bending  moments  and  concluded  that  "the  critical  stress  attained  was  very  close 
to  the  one  of  a  circular  cylinder  subject  to  axial  loading."  (In  Table  2,  the 
critical  stress  and  moment  are  identified  as  «CR>  MCR* )  Consequently, 
although  the  experimental  results  for  models  ICl  and  IC3  pertain  to  axial 
compression  and  "failure  loads,"  an  assessment  is  in  order  to  evaluate  our  own 
modeling  strategy  in  this  light  with  the  objective  of  further  validation  of  the 
employed  modeling.  Very  specifically,  the  experimental  failure  stresses  should 
be  less  than  the  ones  predicted  by  ABAQUS.  How  much  smaller  they  should  be 
depends  on  the  exact  position  reached  (axial  compressive  strain  at  failure)  at 
the  end  of  the  test.  Furthermore,  the  stresses  and  bending  moments  obtained  by 
the  finite  element  procedure  should  not  exceed  the  values  by  bifurcation  theory 
for  purely  elastic  cases  (within  round  off  errors).  In  the  case  of  plasticity, 
the  bending  moment  obtained  by  simple  beam  theory  offers  such  a  guide.  A  final 
point  in  this  discussion  concerns  the  fact  that  the  experimentally  obtained 
stresses  represent  "failure  stresses"  and  not  stresses  at  the  exact  point  of 
instability,  hence  the  drop  in  the  load  that  could  be  carried  there. 

Finally,  we  observe  from  Table  5  that  the  experimental  end  panel  failure 
stress  for  IC3  for  pure  compression  is  36.404  ksi.  Table  6  gives  the 
theoretical  buckling  stress  for  model  IC3  as  115.928  ksi.  At  this  point, 
observe  that  the  differences  between  both  experimental  failure  stress  (25.671 
ksi  and  36.404  ksi  for  models  ICl  and  IC3)  and  the  computational  results 
(49.18  ksi  and  51.72  ksi,  respectively)  are  due  to  the  highly  unstable 
postbuckling  behavior  of  the  structure. 
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There  are  small  differences  between  some  of  the  constants  of  Table  6  and 
their  values  on  Figure  18  (model  IC1)  and  Figure  35  (model  IC3).  This  is 
because  Table  6  was  manually  generated. 

In  conclusion,  in  spite  of  the  differences  in  stress  levels,  which  might  be 
attributed  to  the  fact  that  some  of  the  material  constants  were  inconsistent 
(see  Ref.  26),  it  has  been  demonstrated  that  plots  of  the  critical  parameter  for 
short  axial  wavelength  (local)  buckling  do  indeed  capture  the  mechanism  as  it 
grows  with  load  increase.  The  hoop  stress  distribution  changes  shape  radically 
(compare  Figure  15  with  Figure  22).  Furthermore,  the  cross-section  flattening 
parameter  plot  shows  a  discontinuity  at  the  point  of  short  wavelength  buckling 
initiation. 

SUMMARY 


The  present  strategy  is  easy  to  grasp  and  it  can  be  used  in  conjunction 
with  a  general  nonlinear  finite  element  computer  program  to  predict  the  mode  of 
failure  of  a  cylindrical  shell  subject  to  end-bending  moments.  However,  a  number  of 
problems  exist  that  must  be  addressed,  such  as  verification  of  the 
moment-curvature  relations  for  ring-stiffened  shells. 
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FIGURE  2.  (A)  CYLINDRICAL  SHELL  (SHOWN  AS  A  REAM)  OF  LENGTH  L  SUBJECT  TO  SHEARING 
LOADS  F.  PRODUCING  A  BENDING  MOMENT  M  -  Fa  (THIS  CONFIGURATION  APPLIES 
TO  MOOELS  10A,  ISA.  20A.) 

(B)  CYLINDRICAL  SHELL  OF  A  REDUCED  LENGTH  LR  SUBJECT  TO  AN  EXTERNAL 
BENDING  MOMENT  M  -  Fa  (THIS  APPLIES  TO  MODELS  10AM,  ISAM  AND  20AN, 
RESPECTIVELY.) 


MOOELS  10A  AND  10AN 


MEAN 

CURVATURE 

(1/FT) 

0.000000002700 
0.27777—22-02 
0.393399412-02 
0.92—71712-02 
1.7NUMX-02 
0.— 3594792-02 
0 .073019322-02 
0  •  9444—102-02 
0 .1019— 772-01 
0.109339902-01 
0.1132—442-01 
0.122229532-01 
0.129170732-01 
0 . 1  Ml  19492-01 
0 .142000702*01 
0.1—003742-01 
0.170209292-01 
0 .107920990-01 
0.197249242-01 
0 .202207992-01 
0.20—49992-01 
0.211449712-01 


0.000000002700 
0. 99939—92 -02 
0.— 71912-02 
0.709399792-02 
0.919—1932-02 
0.112902372-01 
0.139937992-01 
0.134179992-01 
0.173009742-01 


BENDING 

MOMENT 

(KIF-FT) 

0.000000002700 
0.393777472702 
0.101—9992703 
0.10742—12703 
0.109721292703 
0.111221742703 
0.112347—2703 
0.1193—492709 
0.11—94942703 
0.11479992C709 
0.113—9712709 
0.113929992703 
0.11—32—2703 
0. 11994—— 70S 
0.11721—72709 
0.11737—22703 
0.119997—2703 
0.11—49292701 
0.11—2—92701 
0.1201143—701 
0.120  —0722701 
0.120492942701 


0.000000002700 

0.71—47422702 

0.107399192703 

0.109—3402703 

0.11—39792701 

0.113100492703 

0.119—1912701 

0.117—7992701 

0.11—24172703 


MOOELIOA 


MODEL  10AN 
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MODELS  ISA  AND  16AN 


MEAN 

CURVATURE 

(1/FT) 


BENDING 

MOMENT 

(KIF-FT) 


I. 

I  UIMMII-H 
•  IMiTUH-H 

a.  ma7taM-ai 
a«uii7tN-a< 
a.M7»4atE-u 
a  >aMii*7c-M 
anai>au(-<i 
a  u«st4(M-ai 

aiiaa  7at7i-aa 
a  c*M7*7«c-ai 
a  iaiia»M-M 
a  ftainui-ai 
a.  *i«iasaac-M 

a  tMMUM-U 
a.  MMN7K>U 
a.  Miuiiw-u 

a  47(4117*8-81 
a  M«MatH-H 
a.  7a*47«*n-aa 
a  7*47*1 t4i-*i 
a. 73*747771-1* 
a.  7t«auu(-u 
a.7tfaM)K-ai 
ammtti-M 
a  7tfni4t(-ai 
a  ai4tit«MHi 
a  a)H4i)ti-u 
a.  wian>*7c-aa 
a.  a«ia74t»-a< 


a.iaaiittii4ai 
a. i347iiaac+aj 
a.  u«u(m»ii 

aiH<u4X7a] 
a  i4tH4aa(Hi 
a<4*i4*7X4ai 
a.  **787ai4i48i 
a  itii444Uiai 
a  NiH7a«H) 
a  i*i?4aaac+ai 
aiili714U4«} 
aiU7f  M4Ii«] 
a  K1M44H4I] 
ltSiaMf4I4|) 
aa*4i4ia*C4ai 
a. «t48S4itc+ai 
a.  nui77«t4«i 

a.  n>aaiiac4ai 

a  H411HHH! 

a  n44iii«4a] 
a  **4*7ia«c+B) 
a.  **711*3714*1 
a  **74**331483 
a  **7*17*41481 
a.  i*ai3884i4si 

a  nai4747C4ai 
a.  naaia*«K4a> 
a.Ma7iJNC4ai 
a  naaaiasi+B) 

a  na4474ji4ai 
a  n *84738(487 


MODEL  ISA 


a  KUIHII-M 
8.  341 1  l*t*C-** 

a  laaaaiaac-M 

a.4i444fm-ai 
a.  U77a  law -aa 

a  444478841 -ai 
a  8—174881  81 


a  iiaaaii4c*ai 
a  114418141483 
8.  1384734*1481 
8  131471841483 
8  147111881481 
■  11 3187141481 
8  11 7871 111481 


MODEL 16AN 
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MODELS  20A  AND  20AN 


MEAN  BENDING 

CURVATURE  MOMENT 

(1/FT)  (KIP-FT) 

O.OOOOOOOOE+OO  0.00000000E+00 

0 . 33794293E-0 3  0 .32341548E+02 

0.12734440E-02  0 . 19720113E+03 

0.22129490E-02  0 .33730443E+03 

0 .31 SO 491 3£ -02  0 . 39873380E403 

0.39473738E-02  0 .42O90343E403 

0.41317S7CE-02  0 . 42434280 EVO 3 

0.441C1407E-02  0 . 42743844E40 3 

0.44303242E-02  0 . 430303S4E+03 

0. 488491 12E -02  0 . 432430 26E+0 3  MODEL  20 A 

0.S11929S3E-02  0.4340457X+03 

0.S3S3C799E-02  0 . 43333480E403 

0 .33880 444E-02  0 .43498942E403 

0 . 3922431 3E -02  0 .43834884E403 

0.40349429E-02  0 .43943071E403 

0.4291 2297E-02  0 .44049992E403 

0 .432341 49E -02  0 . 441 7S433E+03 

0.47400049E-02  0 .44282471 E403 

0.49943932E-02  0 .443481 70E403 

0.72297484E-02  0 .44431007E+03 

0 .7337470 9C -02  0 . 44S00220E403 


O.OOOOOOOOE+OO 
0 . 22222274E-02 
0 . 34111393E-02 
0 . 3400 3830 E -02 
0 .402791 97E -02 
0 . 43730 322C-02 
0 . 47222944E-02 
0. 9244224 9E-02 
0 . 93747342E-02 
0 . 99092334E-02 
0.1 0042273E-0 1 
0. 10121 231E-01 
0 . 10200233E-01 
0. 1027921 4E-01 
0 .103391 99E -01 


O.OOOOOOOOE+OO 
0.30878234E403 
0 .37701491E403 
0.37X319X403 
0 . 39479930E+03 
0.3991 473X40  3 
0.39271 204E+0 3 
0.40 39988 7E40 3 
0.40403331E+03 
0.4039924X403 
0.40 393041 E403 
0.4039194X403 
0.40399010E403 
0. 40X374X403 
0.4037970X403 


MODEL  20AN 


FIGURE  S.  COMPARISON  OF  BENDING  MOMENT  VERSUS  MEAN  CURVATURE  FOR 
MOOELS  SOA  AND  JOAN 
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FIGURE  6.  MOMENT  VERSUS  ROTATION  FOR  MODEL  10AN 


TOTAL  MOMENTS  ARE  FOR  180  OF  CYLINDER 
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NODE  NO 


NSWCTR  86-326 

TOTAL  MOMENTS  ARE  FOR  180°  OF  CYLINDER 


STEi-  6 


SI  £P  4  site  5 

FORCE  2  2  ARM  XX  FORCE  17  ARM  X*  FORCE  17  ARM  XX 

301  0.  30l479E-*04  0.  592331E+01  0.  304174E+04  0  o02250E+01  O.  306793E+04  0  6l204bE«-01 

302  0  120650E«05  O  5879BBE+01  0  121733E+05  0  597946E+01  O  1227BBE+05  O  6G7786E+01 

303  0  6029 11 £+04  0  575037E+01  O  608350E+04  0  585105E-KJ1  0.  613642E+04  0  595065E+01 

304  0  1202S9E+05  0.  553624E+01  0  121307E+0S  0  563861E+01  O  122353E+05  0  574004E+01 

305  0  600720£*04  0.  524067E+01  0  60S581E-f04  0  534508EMH  0  610830E+04  0  544874E+01 

306  O  1 19596E+05  0.  48667 1E+01  0  120334E+05  0  497330E+01  O  12122QE+05  0  507932E»01 

307  0  594796E+04  0.  442240E+01  0  597627E+04  0  453105E+01  0  600905E+04  0  463932E*01 

308  0  118314E«05  0.  391296E»01  0.  U8743E+05  0  402338E-*0l  0  119157E+05  0  413357E»01 

309  O,  593475E+04  0  334692E+01  0  593568E+04  O  345870E+01  0  592712E+04  O  357035E*01 

310  0  111729E+05  0.  273405E+01  0  114054E*05  0284673E+01  0  115142E*05  0  295935E+01 

311  0  490622E+04  020B552E<01  053B661E+04  0219870E+01  0580640E+04  023tl85E*0l 

312  0  455724E-*04  0.  14116BE*01  0  525773E+04  0  15250BE+01  0  594174E+04  O  163B45E+01 

313  -O  783798E+02  0.  7259668*00  -0  962B02E+02  0  839412E+00  -0  125566E+03  0  952836E*00 

314  -0.  48523BE*04  0.  390707E-01  -0  561892E+04  0  153300E»00  -0.  634439E+04  0  266707E+00 

315  -O.  500499E+04  -0.  6350B5E+00  -0  549745E*04  -0  S21645E-*00  -0  591537E+04  -0.  40B228E+00 

316  -0  U2339E+05  -0  128544E+01  -0  114537E+05  -0  U7188E*01  -0  115742E+05  -0  105B30E*01 

317  -0.  593465E+04  -0.  190080E+01  -0  593223E+04  -0  178684E+01  -0.  592279E+04  -0  167279E*01 

318  -0  11SG20E+05  -O.  246984E+01  -0.  1 1B405E+05  -0  235509E+01  -0  1 1B780E+05  -0  22401 6E+01 

319  -O  591970E+04  -0  298253E*01  -0  594356E+04  -0  286663E+01  -0  597292E+04  -O.  275040E»0l 

320  -O  118806E+0S  -O  34301 2E *01  -O  119429E+05  -0  331276E+01  -O  120210E+05  -O.  319488E+01 

321  -O.  595653E+04  -O  38071 5E+01  -0.  599837E+04  -0  36S818£f01  -0  604302E+04  -O  356849E+01 

322  -0  11911 1E+05  -0.  410535E*01  -0  1199B9E+05  -O  398481E+01  -0.  120858E*05  -0.  3B6336E+01 

323  -O  596S19E-»04  -O.  432149E+01  -0  600948E+04  -O  419960E*01  -O  605251E+04  -0  407671E+01 

324  -O  119321E*05  -0.  445227E+01  -O  120200E+05  -0  432953E+01  -0  121051E+05  -O.  420563E+01 

325  -0.  298056E+04  -0.  449613E+01  -0.  300220E+04  -0  43730BE+01  -0  302308E+04  -0.  424B83E+01 


TOTAL  OUT  OF 
BALANCE  FORCE  21 

0  153809E-01 

0  610352E-02 

0.  371094E-01 

TOTAL  MOMENT 

0  699495E+06 

0  705810E+06 

0  71 1425E+06 

FORCE  22 
ON  TOP  QUARTER 

0  954521E+05 

0.  973544E+05 

0  990638E+05 

FORCE  22 
ON  BOTTOM  QUARTER 

-0  954521E+05 

-0  973544E+05 

-0  990638E+05 

POINT  OF  2ER0  FORCE 

0  900390E+00 

0  101B20E+01 

0  1 13577E+01 

PERCENTAGE 
OUT  OF  BALANCE 
FORCE  22 

0  00 

0  00 

0  00 

mi  i««)i  maun 


NOTE:  MODEL  10AN  WITH  IMPERFECTIONS 

STEP  1,  INCREMENT  5,  WITH  PLASTIC  DISSIPATION  -  6686.6  Ib-in,  ANGLE  OF  ROTATION  -  1.217  deg 
TOTAL  BENDING  MOMENT  -  1.30686x10®  Ib-in 


FIGURE  7.  (CONTINUED) 


NSWC  TR  86-326 


sn.e  no 


6  FOR  STRESbt  S 


ANCLE 

LONGITUDINAL 
MEMBRANE  STRESS 

(DEGREES) 

(PSI  ) 

0.  000 

0  627772E+05 

7.  576 

0  623431E+05 

15.  151 

0  61 9201E+05 

22.  704 

0  605740E+05 

30  241 

0  592030E+05 

3/.  761 

0  56B4B2E+05 

45  246 

0.  544551 E+05 

52.  707 

0  5 1 4739E+05 

60  159 

0  405925E+05 

6  7  587 

0  452500E+05 

75  000 

0  453430E  »05 

82  414 

0  207462E+05 

89  817 

-  7549blC+G4 

97  222 

-  354746E  +05 

104  642 

-.  616167E  +  05 

1  12.  067 

-  606059E+05 

1  19  514 

608H0/E+05 

126  990 

-  59 1  Of) /E+05 

134  482 

-.  5/5/9VL  *05 

142  003 

-  557752L+U5 

149  564 

-  53955 /L +05 

157  146 

-  524  W*4L  « 05 

1 64  747 

-  50094 bE  *05 

172  373 

-  502908E  *05 

1  BO  000 

-  49725l»r»u5 

NOND I  MENS 1 QNAL  MEMBRANE 
longitudinal  STRESS 


HOOP  STRESS 
(PSI  ) 


NOND I MENS I ONAL 
HOOP  STRESS 


0  808620E-01 
0  B0302BE-01 
0  797579E-0 1 
0  70024  IE-01 
0  7625B2E-0 1 
0  732250E-01 
0  701424E01 
O  663024E-01 
O  62591 0E-01 
O  502855E-01 
0  584053E-01 
0  267227E-01 

-  9724 7 5F -02 

-  45694  IE-  01 

-  79  1672E-01 
780652E-01 

-  /04192E-01 

-  /61367E  01 

-  /41674E-01 
7 1 84P9F-0 1 

-  6949921.-01 

-  675190E-01 
655566F  -01 

-  647785E-01 

-  64050BE-01 


0  133165E+05 
0  1 226G2E+05 
0  11  19561:  *-05 
0  825290E+04 
0  526420E+04 
0  104852L+04 

-  31907 1E+04 

-  7659B1E+04 
1213B7E+05 

-  155863E+05 
1 92195E+05 

-  21002  1  E+05 

-  223950L+05 

-  212457L+05 

-  20 1 1 64E+05 

-  162286E+05 

-  1231 CUE +05 

-  /65263E+04 

-  3023551 +04 
0  110424E+04 
0  537192E+04 
0  B2L540C  *-04 
0  1 1 1 154E+05 
0  121446E+05 
0  1313651+05 


0  171 527E-01 
0  1 56025E-01 
0  1442G9E-01 
0  106J05E-01 
0  67008 IE-02 
0  135058E-02 

-  410989E-02 
-. 986644E-02 

-  156H56E-01 

-  200763E-01 

-  247563E-01 
270524F -0 1 

-  288465E-01 
273661E-01 

-  259115E-01 

-  209038E-01 

-  15G573F.-01 

-  985719E-02 

-  389457E-02 
0  1 52539E-02 
0  691946E-02 
0  106337E-01 
0  143175E-01 
0  156432E-01 
O  169209E-01 


l  OWCI  f  uCIMHi  *T»r*S  4MC4.LAR  iKttlON 


mo 5*  .•■r.. 


MODEL  10AN 

STEP  1,  INCREMENT  30,  PLASTIC  DISSIPATION  =  46638.8  lb-in,  ANGLE  OF  ROTATION  =  3.008  deg 
TOTAL  BENDING  MOMENT  -  1.42349x10®  lb-in 

FIGURE  8.  AXIAL  AND  HOOP  STRESS  DISTRIBUTIONS  AT  MIDLENGTH  FOR  MODEL  10AN 


INCREMENT  NO.  STRAIN  ENERGY  EXTERNAL  WORK  DONE  PLASTIC  DISSIPATION  ABAOUS  STEP  ABAOUS  INCREMENT 


NSWC  TR  86*326 
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FIGURE  10.  MOMENT  VERSUS  ANGLE  OF  ROTATION  FOR  MODEL  20AN 
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FIGURE  12.  MAGNITUDE  OF  AXIAL  HARMONICS  (OF  COSINE-COSINE  TERMS)  FOR  MODEL  20A 
WITH  IMPERFECTIONS 
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FIGURE  14  DEFORMED  AND  UNDEFORMED  PROFILES  FOR  MODEL  IC1  AT  STEP  3.  INCREMENT  30 
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FIGURE  15.  (A)  DISTRIBUTION  OF  AXIAL  STRESS  AT  MIDLENGTH  FOR  MODEL  IC1 
(B)  DISTRIBUTION  OF  HOOP  STRESS  AT  MIDLENGTH  FOR  MODEL  IC1 
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FIGURE  17.  BENDING  MOMENT  COMPUTATION  THROUGH  CONSISTENT  FORCES  AT  MIDDLE  FOR 
MOOEL  IC1  (FOR  HALF  ONLY) 


FORCE  II  ARM  XX  FORCE  II  ARM  XX  FORCE  II  ANN  II 


lliilll! 

0000000000000^(j)^(jMj)^(jMj)(jMj) 

m*mm*$m$m*m** 

OOOOOOOOOOOOOQ(j)Cj)ij)(j)(j)(j)(j)(j>(ji(j)(j) 

!!!!!««!???«!««!!!!! 


8*888|||H8888 

r>Q+~Q-QC*ri~J>~ 


ooooooooooooooooooooooooo 
i  i  »  i  i  i  i  i  i  i  T  i 


?oo???o???5oo53o3??ooo?oo 

8888  8#  imagine  i8MK|81!|N!!!ie 

8-tf)rv*>r'4aniNil>4«nnni9n#«#3n'«** 
Qfs«fs^rt5^mnpni»opf3Qnici»n»«Q- 
-ft*9‘*N«-n*>Q4frr>gAftNKgir'nNrtgfr 
«n-o~*>®nnini*rv^4(S«n(vna*i~4nn 

oooooooooooooo<po<j>oo<^><pcp<j>(j>^> 

8888888888888888888888388 

5""£*&*F-**23S*?"*8£S&8S 

8>rtr»rtin<«frfvn8>irinnr<rv(hqa^^oi5Pa3S 


_  ..  -  .  _  _*  n  n  8>  4  c*#  n  «r*  n  *  a  •  l  ._  _  „  _ 

r<f»Snininir8---*Na-nrv-.-wf5nio»cM?5^ri( 

ooooooooooooooooooooooooo 

I  I  I  I  I  I  I  I  I  I  I  I 


?  ?  !  !  f 

y  y  yl  in  HI 

S!  (  ;< 

no  n  n  n 


SSooooo3oooo8S3SSoo333ooS 
♦  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

? 

N 

O 

♦ 

O 

o 

♦ 

O 

o 

♦ 

o 

1 

88883*888*8888  88  88*888881? 

-4<Mfrn£Ci~o-*--*»nr»n^*g?>c8a*fr-* 
rv^ofn-coqo^®0‘OQi-o>  4ooor»-*®N  4 
—  -*#<«rNQf#naonQfv^2®-.  —  «  ♦ 

Ui 

N 

5 

8 

♦ 

8 

* 

0* 

-0 

8 

♦ 

•0 

.  1  l  1  11  ™  l  1  ™  1  1  Ml  H  11  "  11  41  MJ  1  >  ~  »  J  ^  M  II  U  u 

r^a^a^rAnmmn  —  —  a  —  — r’lnmnrv^a^frrt 

- 

s 

7> 

s 

- 

000000000000^^9^(j)(j)^00(j)ij)(j) 

o 

O 

o 

? 

o 

Sn^Ss8nSnnnnnninnnnnS3n3 

S~ 

»r  a 

z 

t 

^  a 

N  Ui 
►- 
ui  ac 

^  a 

*i  ui 
►» 

o5 

ui 

u 

8 

Ui  Ui 

n 

8b 

? 

S3 

si 

o 

2  ^ 

o  a 

»-  IN 

i 

*S 

"s 

ui 

N 

88 

i 


IvKWWW 


FIGURE  17.  (CONTINUED) 
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M01  *  VdTSSJfl  H  rv'H  M^IOM  IWMM^iy.q  A0H3N^  MIVfcMS  10 1  '  lOW/U 


UU«M!BMillWWWlMrai«UM*»U» 


NSWC  TR  86-326 


O  O  C  o  -  i‘«  »\ 

o  o  o  o  o  o  c 

4  *  4  4  *  4  ■* 
UJ  Ui  Ul  UJ  UJ  UJ  in 
O  O  O  o  O  -0  V 
ooooNinN 
o  o  o  o  m  ^  h 
oooo-con 
ooocino- 
O  O  O  O  CM  «  - 


ooooooooo 


r-  >n  • 

O  -«  V 
CD  -o  « 


Uj  UJ 
T  ST 

rv  o 

V 

r-  o 


•M  O 

o  o 


II 

£ 

V 

z 

UJ 

n 


n 

X 
o  - 
z 

O 

z 

UJ 


0  0  9  0  0 


o  >n 
ru  -o 
n  n 


o 

z 


o  c  o  o  o  o  o  o  o  o  o  o 


4  4  ♦  4  4  4  4  . 
uiliiUJUJUJUJUJUJ 

no^a5-.n«f«T 
«f  O  r.  -C  o  D  4  £► 
^  f u  m‘«  CD  n  m  n 
■ «  :j  n  Off  t  ^ 
i*3  in  u")  in  a")  ul  u*>  if> 

OODCOOCO 


•\  r u  u  m  oa  cj  .-j  fj 

oooooooo 
♦  ♦♦♦♦»♦♦ 
uj  uj  uj  uj  ui  uj  u;  u 
•o^c^ncDj-oo 
nCDOO—  CDi>f^ 
^  ^  ?  s  Cj  £  r.  - 

..  o  .>  y  o  o  —  r«j 

■j  .*n  n  '"i  cj  .n  o  ^ 

00000900 


o  o  o ooooooooo 

♦♦4*4  +  44444 
JjLuluO'uiUJUJiJU.'UJUJUI 

omocii-niM^MnoN 
o-oo-iODDai-ot  «mn 

^  o  r:  N  O  *f  CM  o  OD  -fl  y 

iiMiin»ooi>>cDiim) 
/,:*.*/>«■  O  O  O  O  O  O  O 
—  rj  — 

ooo ooooooooo 


SE  - 

--  U. 
K  S 
<  - 
>  - 
x 
o 
o 

z 

< 


c 

z 


!???$?*$$??* 

r^mcncprjnjNo^NUjCD 

NftlftflONtWCDtHffl 

rf>pcn^torsNrNNrvrvr» 
cun  - 

666000606060 


n^^^mminininin^r) 
00  0000000000 
♦  ♦  ♦♦♦♦♦♦♦  +  ♦♦ 
UJ  UJ  UJUJUJUIUJliJUJUJUJUJ 

frCD«*m~inn*i>n0‘0* 
CMin^u-^h-N^^-cthin 
-*ra«-nji>-ono>o-<r'-sr 
cm  -*  »nn 

or>  tsoaiH««HHHH 
ncMnjm^njnirucvirarjrj 

00  6000000000 


_J 

c«  - 

-*000 

0 

0 

0 

0 

0 

O 

< 

0  0 

0000 

0 

0 

0 

0 

0 

O 

0 

1  i 

1  4  ♦  ♦ 

+ 

♦ 

+ 

4 

4 

4 

UJ  UJ 

UJ  UJ  UJ  UJ 

UJ 

UJ 

oJ 

UJ 

UJ 

UJ 

K 

in  43 

n  m  0  0 

«• 

Is- 

in 

*• 

CD 

in 

i>  0 

r  4  *  0  cm 

CD 

n 

<■ 

•0 

X 

n  0 

sr  n  0  0 

O 

0 

N 

m 

CM 

u 

<  Is- 

0  m  in  Is* 

-0 

in 

n 

CM 

0 

r 

n  n- 

N  if)  N  J) 

0 

UJ 

>0 

-0 

<3 

-0 

v 

0  N 

cd  -  m  -0 

•0 

-0 

>0 

<0 

O 

-0 

0  c 

0000 

0 

0 

0 

0 

O 

0 

-  c 
o  o 
1  4 

UJ  UJ 
0*  f  0 
O  >0 

n  -© 

n  rx 
in  ru 


o  o 
$? 

§2 
o  u 
o  m 

iD 

-  n 


00000000 

44444444 

UJUJLUUJUJUJUJUJ 

-nthriiMDON 

inconj-o^r  —  ny*) 
r'tf'>srrootD-o<r 
«r  Ch  Ch  J  (md  a)  a) 
*0000000 


OO  OOOOOOOOOO 


—  OOO 
OOOO 


O  O  O  t 

O  O  o  > 


44444444 


*  CO 
l>  - 

n  cd 

*>  Is 


-*»  CD  * 
«  CD  ^  CD 

-«  —  n  m 
•n  -o  in  43 
Is  m  o  o 
ca  ~  1 n  N 


Ui  UJ  UJ  UJ 
03  30  -0  (h 
-OO  N 
o  'O  m  n 
n  o  n  * 
0  o  a  a* 
n  r*.  >c  -o 


UJ  UJ 
•C  CD 
nj  -• 
n*  ro 
-  0* 
o  0 
*0  -0 


00  OOOOOOOOOO 


*  *  */  n  r  ii  n  n  n  n  n  n 

00  ocoooooooo 

I  .  I  I  I  I  I  I  I  I 

UJ  CL  xaJUJUJUJUJUJUJUJUJ 
nj  ®  TNncprjp-anr) 
nm  o**CMOfM*rM«c.o 
in»n  -.*cDn(>ir®f6-*{M 
©  ftf  -‘ON’OttlMftND'O 

00  omflo^^ooooo* 

—  0*  ^^-CCDOCDaDCDrsrv 
00  OOOOOOOOOO 


-  I'J  O  ♦  IT)  -0  N  B  0  O  -  fil 


3 

z 


o 

o 


ec 

D 

o 


30 


CROSS-SECTION  FLATTENING  PARAMETER 
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FLATTENING  PARAMETER  -  MEAN  CURVATURE 


FIGURE  19.  CROSS-SECTION  FLATTENING  PARAMETER  Aw/h  VERSUS  MEAN 
CURVATURE  FOR  MODEL  IC1 
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EXTERNAL  WORK  -  BENDING  MOMENT/ULTIM  MOMENT 


FIGURE  20.  M/MULT|MATE  VERSUS  WORK-DONE  FOR  MODEL  IC1 
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FIGURE  24.  DEFORMED  AND  UNDEFORMED  PROFILES  FOR  MODEL  IC1  (WITH  ADDITIONAL 
CONSTRAINT  W  =  0  AT  END  PLATE)  FOR  STEP  3,  INCREMENT  20 
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CRITICAL  PARAMETER  FOR  SHORT  AXIAL  WAVE  BUCKLINC 
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FIGURE  32.  (A)  DISTRIBUTION  OF  AXIAL  STRESS  AT  MIDLENGTH  FOR  MODEL  IC3 
(B)  DISTRIBUTION  OF  HOOP  STRESS  AT  MIDLENGTH  FOR  MODEL  IC3 
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FIGURE  34.  (CONTINUED) 
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(7)  MOMENT  OF  INERTIA  OF  UNDEFORMED  CROSS-SECTION 
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NOMENCLATURE 


E,F,G,H 


EiJ»  Fi.J* 
GIJ*  HIJ 


F0 

f*(tn,6y  ) 
F(x,y,z) 


CR 


l 

2L 

M 


M 


BR 


Mq 


Matrices  whose  coefficients  are  given  by 
EiJ »  FiJ»  GiJ«  HiJ 


Form  of  coefficients  defined  in  Equations 


and 

respectively. 


(A-41)  -  (A-49) 

(A-50)  -  (A-54) 

( A-55)  -  (A-59) 

(A-60)  -  (A-64) 


Axial  Force  based  on  yield  stress  (2Rhoy) 

Defined  in  equations  (A-30),  (A-31)  and  (A-32) 

Symbolic  representation  of  the  surface  of  the  cylinder 
(F  =  X2  +  y2  -  r2) 

Shell  thickness 

Moment  of  inertia  of  undeformed  cross-section  (ttR-^i) 
Critical  curvature  at  bifurcation  /  E  \  h 


yj  3(  L-v  - )  /  R 


Total  length  of  shell 
Total  Length  of  shell 


Upper  bound  in  axial  wave  numbers  (If  only  a  sinusoidal 
distribution  is  used,  then  there  will  be  M  points,  otherwise, 
M+l) 


Critical  Brazier  moment 


9  yj(  l-v?-) 

Critical  moment  at  bifurcation 


ERh 


Rir 


Plastic  Moment  based  on  yield  stress  (4R“hOy),  where  R  is 
mean  radius 
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N  ”  Upper  bound  In  peripheral  wave  numbers.  (There  Is  a  total  of 

N+l  points  peripherally) 

_n  *  Outward  normal  to  cylindrical  surface 

(Nx,Ny,Nz)  =  Components  of  ^  in  global  Cartesian  system  (x,  y,  z ) 

R  ■  Actual  radius  of  cylindrical  shell  including  imperfections  (In 

formulas  for  oCR»  MCR»  KCR»  MBR»  M0 *  F0»  and  R  ls 
the  mean  radius  as  given  by  Rg) 

Rq  “  Mean  radius  of  circular  cylindrical  shell  (perfect  shell  or 

radius  of  imperfect  shell  prior  to  introduction  of  imperfections) 

(R,0,z)  *  Cylindrical  frame  of  reference 


(x,y,z) 


-Jl-  k 
(M-l)  * 

Global  Cartesian  system  of  reference 

Critical  stress  at  bifurcation  h 

R 


Gradient  of  function  F,  l.e.  , 

(Fxj  Fy ,  Fz)  or 

3F,  jjF, 

3x  3y  3z  / 

Imperfection  coefficients  (usually  of  small  amplitude  compared 
to  R()) 

2tt  u  ,  where  K  is  the  total  number  of  points  Including 
( K-l)  the  0  and  2r  points  (X  is  odd) 

Yield  stress 

Magnitude  of  gradient  of  F,  i.e.  , 

?  ?  7  i 

(F  '  +  e  +  F  Y* 
x  y  z 
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APPENDIX  A 


ANALYSIS  OF  A  NEARLY  CYLINDRICAL  SURFACE 


In  this  section  we  consider  the  definition  of  a  nearly  cylindrical  surface 
based  on  a  cylindrical  frame  of  reference  (R,9,z).  The  object  Is  to  obtain  the 
perturbed  cylindrical  surface  from  an  ideal  one  of  mean  radius  R0  and  total 
length  2L. 

Two  types  of  perturbations  will  be  considered  for  the  imperfections: 

(1)  No  imperfections  at  z  c  0  and  z  ~  21.. 

Symmetrical  imperfections  peripherally.  Therefore,  the  local  radius  R 
is  given  by 


M  N 


m=l  k=0 


where  <5mlc  are  the  imperfection  coefficients.  The  equation  of  the  surface  F  of 
the  thin  hollow  cylindrical  shell  is: 


F(x,y,z)  =  +  y-  -  R2  = 


M  N 


and 

x  =  RcosO 
y  =  RsinQ 
z  =  z 


( A-2 ) 


(A-3) 


Next,  we  need  to  define  the  outward  normals  n,  =  (nx,  ny ,  nz)  to  the 
perturbed  surface.  This  can  be  accomplished  by  obtaining  the  gradient  of  the 
function,  which  defines  F(x,y,z). 


NSWC  TK  8o-32b 


We  must  obtain  7Ff  wh .?rK 

-  -(i  •  f .  ri-G.-vO 

IVFI  -  <■  r  ►  F‘)  l/2 

‘  '  (’*  '  ^  '  ’*  )*(F«  '  %  •  OA 


(  4-4 ) 


t  \-5; 


( 4-6 ; 


To  calculate  the  components  of  VF  we  employ  the  form  of  the  function 
given  by  Equation  (A-2),  i.e. 


f  -  2x  -  2:3 

oX  3x 


2y  -  2K— 

ax 


^  =  -  2:3 


(4-7; 
( 4-8) 
<  4-9) 


— 

R(  z , 

0)  and  : 

since  z 

Is  i 

nde| 

3R 

3z 

.  JR 

30 

30 

a  r 

~ 

3/. 

3x 

30 

3x 

3x 

JR 

+  3* 

30 

3  0 

JR 

3z 

ay 

an 

ay 

ay 

3  0 

vA-T.i; 


( A— 1 1  ) 


But  9  =  cot 


-[?] 


( A- 1 2 ; 


( A  - 1  3 ) 


3  0  m  cos  0 
3y  R 


( A- 1 4 ) 


Substituting  Equations  (A-13)  in  (A-10)  and  (A-14)  in  (A-ll)  we  get 


3R 

sinfi 

3  R 

3x 

R 

30 

3R 

COB  0 

>R 

,3y 

R 

3  0 

Equations 

(A- 7 ) 

- 

f  -  2R 
3x 

£c0  8  0 

+ 

( A  - 1  5  ) 


( A- 1 6 ) 


t-9)  become,  in  view  of  (A-15),  (4-16)  and  (A-3). 
s  i  n  n  3  R  1  ,  . 


sin  n  3R  1 

R  3  0  J 


(4-17) 


CV.v’A’.VA\V.vWt’«l%\V 


tMHomnma 
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jjF  _ 

2R  £sin0  + 

cos  0 

3y 

R 

3F 

-  2R 

3  z 

From  A. 1 

M  N 

(A-18) 

(A-19) 


3R 

30 


k  6 


mk 


m*l  k*0 

M  N 


Sin[ll  2] 


cos  [k0] 


(A-20) 


_3R 

3z 


(s)  fmk  COS 

mr 

2L  2 

\  / 

_  ... 

m*l  k*0 

and  finally  Equations  (A-17)  through  (A-19)  become 

M  N 

cos  0  -  sin  O 


cos  [ k0 ] 


(A-2 1 ) 


~  -  2R 
3x 


3F 

~  *  2R 
3y 


m=*l  k=0 

M  N 


sin 


mu 

2L 


z  sin  (ko) 


(A-22) 


slnO  +  cos  0 


m“l  k=0 


sin 


13 

2L 


z  sin  (k‘  ) 


(A-2J) 


M  N 


~  -  -2R 

oZ 


.  2L  /  ‘mk  C°S 


2L 


cos  [k0 


( A-24 ) 


m*l  k*0 


and 


I  VF I  - 


(fi) 


mk 


cos 


ni  tt 


2L 


cos  (k-')l  + 


(A-25) 


lin[li;  z] 


2  T  1  /2 


sin  (ktO 


A-3 
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The  axial  Interval  2L  (or  it  )  is  divided  in  (M-l)  equal  segments  t^  - 

2 

The  peripheral  interval  2n  is  divided  in  K-l  equal  angles  0  *  ,  g 

P  (.  K.—  1 ) 

where  K  is  total  number  of  points  including  the  0  and  2n  ends  (K-odd) 
f  for  the  portion  of  Equation  A-l  to  be  represented  by  the  series  and 


C  "  2LZ 


M  N 


5»k  Bln 


[mu  1 


cos  [k.0  ,,) 


m-l  k-0 


M  N 


Z  Z  6”k  Sin  m  Ck  COS  [k°u 


m-l  k-0 
The  coefficients  fm>x  are  given  by 

if  _  i 

n-M-2 


X*  (M-l)  (K-l) 


S  *(tn,euj  sin  cos  ^  • 


n=*l  u'O 


and 


A  .  ,  ,  ,  C<-1) 


m  -  1,2, ,  M-2 


K-l 

n=M-2  u  =  — 


i,o  T*PT)  Tk^H 


S  sln  T^n'1 


n-1  P-0 


(M-l)k 

» 

Setting 


(A-26) 


|^,](A-27> 


(A-28) 
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M-l 


K-l 


6m,  JW  '  (ThH  Tk3!!  ^  22  1  ( 

2  n= 1  .=0 


1  n  ’ 


-i)‘ 


(A-29) 


for  m  •  1,2,...,  M-I 


and 


f  (Cn  ■  "o)  ’  2  f(tn  ’  ll) 

'  (*"  '  Vi)  '  I  '('n  '  '  Kri) 


for  n  t  0  and  M-l 
for  a  f  0,  M  -  1 


(A- 30) 

( A-31 ) 


due  to  the  effect  of  the  trapezoidal  rule  during  integration.  Also  note  that, 
if  the  terms  corresponding  to  tQ  were  needed,  then 


'*('<>  •  9o)  ‘  1  f(‘o  •  "„) 


(A-32) 


and  similarly  for  the  other  3  points 


^M-l  ’  0^’  f  (C0  ’  "Kzl'j’  f  ^M-l  ’  ’Jill] 


( 2 )  General  type  of  Imperfection . 


In  the  general  type  of  imperfection,  the  cylindrical  surface  in  the  axial 
(0,  ■  )  and  peripheral  (0,  2  )  intervals  is  given  by 


M  N 


Ro  +  f( 

z,  -)  3 

xo  + 

EE 

1  F.  , 

|  mK 

m=0  k=0 

+  Fmk 

cos  | 

1  2m  z  ) 

1  sink 

+ 

+  Gmk 

sin 

^2m ■■ z  j 

1  s  i  n  k 

+ 

+  H  L 

mk 

sin 

^2m” z  ) 

,  1 

|  sin  ik'  :  1 

whe  re  t- 

*  2L 

>»(--) 


cos  (k")  + 


(A- 33) 


A  -  r> 
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The  total  number  of  equldLstant  points  axially  is  2M  +  1  and  the  total 
number  of  points  peripherally  is  2N  +  1,  Including  end  points. 

Discretizing  and  setting 


for  p  ■  0,1,2,  ....  2M 


(A-34) 


2N  q  *  U  q 


for  q  *  0,1, 


(A-35) 


We  can  abbreviate 


p)  ’  >(*  ?)  ’  % 


(A- 36) 


so  that 


f (%.',)■  fM  ■  E  E 


E  ,  cos  (ray  )  cos  ( !<  ■  )  + 

mi<  p  q 


ra=0  k “0 


+  F  ,  cos(ray  )  sin(k  )  +  G  ,  sin(roy  )  cos(k"  )  +1 
mk  p  q  mk  q  q  f 

+  Hn,^  sin(rayp)  sln(k-q)  » 

The  Fourier  coefficients  Fm^ ,  Gm|<,  and  Hm!<  can  be  obtained  if 

we  multiply  Equation  ( A  — 3  7 )  in  turn  with  a  product  of  either 

cos(  yp)  or  sln(,yp) 


( A  —  3  7  ) 


cos(*0q)  or  3 1  n  ( •<  q  ) 

and  sura  in  both  directions  (double  suras),  while  invoking  certain  orthogonality 
conditions  such  as 


2M-1  0  for  *  a 

E]  cos(.  ^p)  cos(ra«fp)  "  M  for  *  m  M 

p”0  2M  for  .  =111“  M 

2M-1  0  for  ^  ra 

E  8in(..fp)  sln(rayjp)  =  .'!  for  .  =  in  M 

p*0  0  for  ..  =  a  =  M 


(A- 38) 


i A- 39) 


2M-1 

V.  sln(ra<Pp)  cos(m^p)  ”  0  for  all  ,t 

p»0 


( A- 40 ) 


NSWC  TR  86-326 


They  are  as  follows: 


2M-1 

2N-1 

E0,0  *  _L 

.  l 

E 

E 

( A— A 1 J 

2N 

2M 

p=*0 

q-0 

p=*2M-l 

2N-1  q 

eo,n  *  _L 

.  1 

E 

E  (-D  f(p,q) 

( A-42 ) 

2M 

2M 

p~o 

q=0 

2M-1 

2N-1 

em,o  ” 

.  1 

E 

E  (-D  f(p,q) 

( A— 4  3 ) 

2N 

2M 

p=o 

q=0 

em,n  “  _L 

.  1 

2M-1 

E 

2^-l  p+q 

23  f(p,q) 

(A-44) 

2N 

2M 

P=0 

q=0 

2M-1 

2N-1 

E0,k  *  _L 

.  1 

E 

23  f(p,q)  cosCkt') 

(A-43) 

M 

2M 

p-0 

q=0  q 

k  =  1,2,  ...  N-l 

2M- 1 

2N-1 

Etn,0  =  _J_ 

.  1 

E 

23  f(p.q)  cos(m<pu) 

(A-46) 

2M 

M 

P“0 

q-0 

ra  =  1 , 2  ,  ...  M- 1 

2M-1 

2N-1 

Em,k  ’  ^ 

.  1 

E 

23  f(p>q)  cos(ro<iPn)  cos(k') 

(A-47) 

N 

M 

p;0 

q  o 

2M-1 

2N-]  q 

Ero , N  =  1 

.  1 

E 

23  (“D  f(p,q)  cos(nxpD) 

(  4-48  .) 

2N 

M 

p=o 

q  =  0 

for  m  =  1,2,  ...  M- 1 

p=2M- J 

2N-I 

eM , k  ’  1 

.  J_ 

E 

23  (_1)  f(  p,q)  cog(  k  '  n !> 

'  A-41)) 

2M 

t  r 

p'O 

q  =  0 

for  k  =  1,2,  ...  ,  N-l 

r  m ,  i )  K  ^ 

for  u 

=  0,  1  , 

2,  ...  ,  M 

(A- 30 i 

l'm,M  ” 

for  mi 

=  0,  1, 

...  ,  M 

i  \-  5  1  ) 

2M-1 

2N^1 

Ki,k  =  _L 

.  1 

V" 

E  f(P»q)  cos(^„)  sin(ki  „) 

f A- 52) 

N 

M 

p~o 

q-0 

»  =  1 . 2 . M- 1 

k  =  1,2,  ...  ,  N-l 

2M- 1 

2N-1 

hV*  =  _L 

.  _1_ 

E 

E  f  ^  P * M }  sln(k  '  ) 

(  \-3Jj 
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2M-1 

2N-1  p 

FM,k 

s 

1 

.  1 

E 

V  (-1)  f ( p, q )  sin(kQq) 

(A-54) 

2M 

N 

p=0 

q*0 

Ir  a  1  9 

IV.  1  •••  > 

N-l 

G0,k 

m 

0 

for 

k  -  0,1, 

. .  .  ,  N 

(A-55) 

GM,k 

- 

0 

for 

k  =  0,1, 

...  ,  N 

(A-56) 

2M-1 

2N-1 

Gm,k 

- 

1 

.  1 

E 

£  f  ( p,q)  sin(m<pp)  i 

cos( k0  q) 

( A-57 ) 

M 

N 

p=o 

q*0 

f or  m  =  1,2,  ...  , 

M-l 

k  =  1 , 2 . 

N-l 

2M-1 

2N-1 

Gm,n 

= 

1 

.  1 

E 

V  (-1)  f(p,q)  sin(nnpp) 

( A-58) 

2N 

M 

p-0 

q=0 

for  m  =  1,2,  ...  , 

M-l 

2M-1 

2N-1 

Gm,  o 

= 

1 

.  _1_ 

E 

f(p,q)  sin(m^p) 

(A-59) 

2N 

M 

p=o 

q*0 

for  m  =  1,2 . 

M-l 

Hm,o 

» 

0 

for 

m  ■  0,1,2 

...  ,  M 

( A-00) 

^m,N 

= 

0 

for 

m  =*  0,1, 

.  ..  ,  M 

( A-6 1 ) 

H0,  k 

= 

0 

for 

k  -  1,2  . 

..  ,  N-l 

( A-62 ) 

H:i,k 

» 

0 

for 

k  *  1,2  . 

..  ,  N-l 

(A-63) 

2M-1 

2N-1 

1  *ra ,  k 

- 

1 

.  1 

E 

£  f(p,q)  sin(nxpp) 

s i n( k" q  ) 

\  A-b-t ) 

\' 

M 

P‘0 

q*3 

form-  1,2,  ...  , 

M-l 

k  -  1,2,  ...  , 

N-L 

Kquat  ions 

(  4-41 ) 

-(A-49),  ( 

4-50)- (A-54 ) ,  ( A-55 )- (A 

-59) ,  and  (A-bUj- 

(A -64  ) 

for  tin  matrix  of  coefficients  K,  _F ,  (i,  and  H  ire  displayed  below  i  n.i  i  c  a '  !  ny  l.e 
zero  elements  only.  Tiie  si /.e  is  ( ?fM)(\'+l ) ,  wlicre  2M+1  indicates  the  totii 
number  of  axial  points,  wiile  2N+1  tbe  corresponding  peripner.il  number. 

By  Equation  (A-33)  we  obtain  t’iat  the  equation  ot  tue  "almost"  ci>-.  u'.iri. 
cylindrical  surface  lF(x,y,:)  is: 


nswc  TR  8 b — 3 2 r> 


M  M 


rn-0  k~9 


S  H  Pmk  co8(«»p)sin(k0q)l(  £  ]T  :’mk  sl nC%)cos(k9q) 


rn-0  k  — i) 


H  ,  sin(rn<p  )sin(k0  ) 
nk  p  (] 


(A-69) 


where 


x  =  R  cos0 
y  =  R  sin0 


( A-70) 


Furthermore,  Equations  (A-7)  through  (A-9)  or  (A-17)  through  (A-19)  are  still 
valid  in  defining  the  gradient  to  the  surface  F,  i.e., 


— -  =  2R  cos  0  + 

(?X 


—  =  2R  sin  0  - 

ay 


i£  =  2r  * 

3z  1  DO 


sin  B 

_3R 

R 

30 

cos  9 

_3R 

R 

3  0 

(A-7 1 ) 

(A-7 2) 

(A-73) 


M  N 


1-11“  c°s(~rz)  - 


sin(ko)  +  k  F^  cos/ — zj  cos(kO)  - 


m^U  k=0 


-  k  G^  sin  — z j  sin(ko)  +  k  sin/— — zj  cos(ko) 


(A-74) 


M  N 


(-f)  _Emk  sIn(^fz)  c°s(k">  -  ' 

m“0  k*0 

+  cos^— ■—  z^  cos(ko)  +  cos|-— zj  s 


Fmk sin rrzJ sin(kf,)  + 


+  G  cos 
mk  ' 


sin(ko)! 


(A-75) 
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“ 

M  N 

2R 

„  sin0 

COS0  +  — — 

EE- 

K  | 

m=0  k=0 

cos(m<p  )  sin(k0  ) 

p  q 


+ 


+  F  k  cos(m^)  cos(kG^)  -  Gmlr  sin(m<p^)  sin(k0  _)  + 


+  sinCmc^)  cos(k0^) 


mk 


2R 


sin  6  - 


cos  0 


M  N 

II' 

m=0  k=0 


-E  ,  cos(m<p  )  sin(k'3  )  + 
mk  p  q 


(A-76) 


+■  F  ,  cos(m<p  )  cos(k@  )  -  G  ,  sin(m<p  )  sin(k6  )  + 
mk  p  q  mk  p  q 


+  H  ,  sin(ra<p  )  cos(k0  ) 
niK  p  q 


(A-77 ) 


2R 


>1 

N 

R 

"  y  (-t) 

L 

m= 

j  Z-j  \  1  ) 

0  k-0 

,  cosCmcp  )  cos(kB 
m?:  p 

i+ 

M  N 

:EE(f) 

oi— 0  k=0 

-E  sin(ra<p  )  cos(kP  )  -  F  ,  sin(m<p  )  sin(k0  ) 
mk  n  n  int-'  rt 


q  idk 


IDK. 


(A-78) 


-E  ,  cos(m<p  )  si;i(k0  )+F  ,  cos(m(p  )  cos(kO  )• 
mk  p  q  nu  p  q 


j  ,  sin(m<»  )  sin(k0  )+H  ,  sin(raa>  )  cos(kO  ) 
mk  Tp  q  nik  p  q 


M  :j 

II 

m=0  k-0 


km  a 


-E  ,  sin(mu>  )  cos(k"  )-F  ,  sin(m<p  )  sin(k0  )  + 
me;  p  q  mk  p  q 


cosimcp  )  cos(k  )HI  ,  cos(rn<p  )  sin(k0  ) 
p  q  mk  p  q 


(A-71)) 


il-'i-  i  ii.irm.il  to  the  .surface  is  given  by  Equations  (A-5)  an  ! 
-  .  :  1  pint  ions  (A-76)  through  (A-7')). 
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APPENDIX  B 

SUBROUTINE  TO  RUN  ABAQUS  BEND  BUCKLING  ANALYSIS 

c 

C  MPC  SUBROUTINE  FOR  BE ND-BU CKL I NG  PROBLEMS 

C 

SUBROUTINE  UPC ( UE , A , JDOF , N , J TYP E  ,  X  ,  U ) 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

DIMENSION  A(N)  , J  DOF ( N  )  ,  X ( 6  ,  N  )  ,  U  (  6  ,  N  ) 

DIMENSION  B1GAS(3) ,AS(3) ,A3(2) 

C 

c  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

c  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C  XXXXXXXXXXXXXXXXX  NOTES  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

C  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX XXXXXXXXXXXXXXXXX 

c 


c 

c 


c 

r 

xxxx 

UE 

= 

VALUE  OF  TOTAL  DISPLACEMENT  AT  ELIMINATED  D.O 

F  .  X 

V 

c 

c 

xxxx 

A  (  N ) 

- 

ARRAY  OF  DERIVATIVES  Of  CONSTRAINT  FUNCTIONS 

XX 

c 

xxxx 

JDOF  (  N ) 

— 

ARRAY  OF  D.O.F.  IDENTIFIERS  AT  NODES  INVOLVED 

XX 

c 

r 

xxxx 

IN  CONSTRAINTS 

XX 

c 

xxxx 

N 

— 

NUMBER  OF  D.O.F.  INVOLVED  IN  CONSTRAINTS 

XX 

c 

xxxx 

WHEN  JTYPErl,  N  IS  THE  NUMBER  OF  ELEMENTS  FOP 

XX 

c 

xxxx 

WHICH  THE  VERTICAL  RIGID  BODY  MOTION  WILL  BE 

XX 

c 

xxxx 

REMOVED 

XX 

c 

xxxx 

NOTE  THAT  WHEN  J  TYP  E=  1  ,  ALL  OTHER  INTEGERS 

XX 

c 

c 

xxxx 

APPEARING  IN  THE  MPC  CARD  ARE  N 

XX 

c 

xxxx 

JTYPE 

= 

CONSTRAINT  IDENTIFIER  GIVEN  ON  CORRESPONDING 

XX 

c 

r 

xxxx 

•MPC  DATA  CARD  (FIRST  ENTRY) 

XX 

L 

c 

xxxx 

X  (  6  ,  N  ) 

- 

ARR*Y  OF  ORIGINAL  COORDINATES  AT  THE  NODES  IN 

XX 

c 

r 

xxxx 

THE  CONSTRAINT 

XX 

t* 

c 

xxxx 

U  (  6  ,  N  ) 

= 

ARRAY  OF  IOTAL  DISPLACEMENTS  A I  THE  NODES 

XX 

c 

xxxx 

INVOLVED  IN  THE  CONSTRAINT 

XX 

c 


c 

c 


C 

xxxx 

WHEN 

JTYPE 

=  2. 

THEN  THERE 

WILL 

BE  5  NODE  NUMBERS, I. E. 

XX 

c 

xxxx 

SHELL  NODE  TWICE  (UZ 

,UX)  AND 

XX 

c 

r 

xxxx 

BEAM 

NODE  THREE  TIMES  (UZB.UXB,  PHIZ) 

XX 

V 

c 

xxxx 

WHEN 

JTYPE 

-  3 . 

THEN  THERE 

WILL 

BE  ?  NODE  NUMBERS, I. E. 

XX 

c 

r 

xxxx 

THE 

SHELL  NODE 

PHIY 

AND  THE  BEAM  NODE  PHIYB 

XX 

u 

c 

xxxx 

WHEN 

JTYPE 

THEN  THERE 

WILL 

BE  3  NODE  NUMBERS, I. E. 

XX 

c 

xxxx 

THE 

SHELL  NODE 

TWICE 

(PHIX.PHIZ)  AND 

XX 

c 

r 

xxxx 

THE 

BEAM  NODE 

ONCE 

( PHIYB) 

XX 

1/ 

c 

xxxx 

WHEN 

JTYPE 

=  1  . 

THEN  THERE 

WILL 

BE  ALL  THE  NODE  NUMBERS 

XX 

c 

xxxx 

ASSOCIATED  WITH 

1  1  HE 

BEAM  NODE, INCLUDING  THE 

XX 

c 

xxxx 

BEAM 

NODE 

XX 

c 

c 

C  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

c 


B-  1 


oooo  ooo  ooooooo  ooo 
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c 

IF  ( JTYPE . NE . 1 )  GOTO  200 

XXXX  JTYPE  i  1 

NODES  =  N-  1 
NN  =  ( NODES- 1 ) /2 
UE  =  0.0 

DO  10  *1=1, NODES 
JD0FU1)  =  1 
T  =  1.0 

IF  U1.E0.1  .OF.  K1.EQ. NODES)  GOTO  5 
T  r  4.000 

IF  ((  K1/2)«2.NE.K1)  Ti2.000 
5  CONTINUE 
A  (  K  1  )  =  T 

UE  r  UE  *  T»  U(  1  ,K1  ) 

10  CONTINUE 

A  (  N )  =  -6.0»FL0AT(NN) 

UE  =  -UE  -  A(N)»U(1,N)  +  U(1,1) 

J  DOF  C  N  )  =  1 
GOTO  990 
200  CONTINUE 

XXXX  JTYPE  =  2 

XXXX  SET  UP  DIRECTION  COSINES 

XXXX  INITIAL  DIRECTION  COSINES,  SINCE  PHIO  =  0.0, 
XXXX  ARE  1.0000,0.0000  RESPECTIVELY 

CUIGF  =  1.0000 
SflIGF  r  0.0000 

XXXX  DIRECTION  COSINES  OF  ADDITIONAL  ANGLE  PHIY 

CFIZ  =  COS(  U ( 5  ,  N  )  ) 

Sr  I Z  =  S I  N  (  U  (  5  ,  N  )  ) 

IF  (  JTYPE. NE. 2)  GOTO  JOO 

XXXX  JTYPE  =  2 


C8IGF 
SB  IGF 
-A(  1  ) 
-  A  (  2  ) 


CFIZ  -  SBIGF*SFIZ 
CFIZ  ♦  CBIGF»SFIZ 


(U( 1 , 1 ) 

(  U  (  3  ,  1  ) 

U  (  3  ,  N  )  h 
U( 1 ,  1  )  • 


X(  1  ,  1  ) 
X  (  3  ,  1  ) 


U( 1 ,N) 
U(3,N> 


-  X(  1  ,N) )  •  A  C  1  ) 

-  X  (  3  ,  N  )  )  •  A  (  2  ) 


+  X  (  3  ,  N  )  -  X  (  3  .  1  )  ♦  (  J  (  1  ,  N  )  +  X  (  1  ,  N  )  - 
-  X( 1 ,  1  )  )•  A(2)/A(  1  ) 


JD0F(1)  = 
JDOF ( 2 )  = 
J  DOF ( j  )  = 
J  DOF ( 4  )  s 
J  DOF ( 5 )  = 
GOTO  99° 
300  CONTINUE 


N  '  N  A  A* 


W 
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If  (JTYPE.NE.3)  GOTO  400 
C 

C  XXXX  JTYPE  =  3 

C 

A  (  1 )  s  1.000 
A (2)  s-1.000 
JD0F(1)  =  5 
JDOF ( 2 )  =  5 
UE  =  U ( 5 , 2  ) 

GOTO  990 
400  CONTINUE 
C 

C  XXXX  JTTPE  =  4 

C 

A (  1  )  =  CBIGf»CfIZ  -  SBIGfSFIZ 

A( 2 )  s-CBIGF*  SF I Z  -  SBIGF»CFIZ 

A  (  3  )  =  U  (  4 , 1  )  •  A  (  2  )  -  U  (  6  ,  1  )  •  All) 

JDOF ( 1 )  =  4 

JDOF ( 2 )  s  6 

JDOF ( 3 )  =  5 

UE  =  -U ( 6 , 1  )  •  A ( 2 ) / A (  1  ) 

990  CONTINUE 
C 
C 

RETURN 

END 
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